Is the Adiabatic Approximation Inconsistent? 
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In recent Letters Marzlin and Sanders and Tong et al. || study an adiabatically varying Hamil- 
tonian h(t) that generates the time evolution U(t) and its dual H(t) that generates the evolution 
U it). Marzlin and Sanders show that inconsistent results are obtained if an adiabatic approxi- 
mation is used to calculate H(t); Tong et al. show that the adiabatic approximation can be very 
inaccurate when applied to the exact dual Hamiltonian H(t) even if it is an excellent approximation 
for h(t). We show that these two observations are equivalent and are not inconsistent with the 
adiabatic theorem because in general, even if h(t) satisfies the conditions of the adiabatic theorem, 
H{t) will likely violate those conditions. 



The adiabatic theorem is the basis of an approxima- 
tion scheme that was discovered at the dawn of quan- 
tum mechanics Q and that has been in widespread 
and continuous use ever since. Applications range from 
two-level systems (such as nuclei undergoing magnetic 
resonance or atoms interacting resonantly with a laser 
held) to quantum held theory (where a low-energy ef- 
fective theory is derived by integrating out fast, high- 
energy degrees of freedom). Two decades ago, Berry 
uncovered the beautiful geometric structure underly- 
ing the adiabatic approximation , leading to a resur- 
gence of interest in the subject and to new applications 
Q, ■ More recently, it has been proposed that Berry 
phase effects lead to quantum phase transitions that lie 
outside the usual Landau-Ginzburg- Wilson paradigm 
[H. The adiabatic theorem is also the basis of a newly 
proposed quantum computing scheme [(|. Consider- 
ing the significance of the adiabatic approximation to 
quantum physics, the discovery of an inconsistency 
would be most disturbing. In a recent Letter Mar- 
zlin and Sanders Q ask whether such an inconsistency 
might exist, at least for a class of Hamiltonians. That 
question has been further studied by Tong et al. |8| 
and subsequently commented on in Refs [1G3- The 
purpose of this Letter is to show that there is no in- 
consistency. 

Refs [3| and [8j start with a time-dependent Hamil- 
tonian h(t) for which it is presumed that the adiabatic 
approximation is accurate. The evolution operator for 
h{t) is denoted by U(t), the solution to 

i^-U(t) = h(t)U{t), (1) 

with U(0) — I, where I is the identity operator. Next, 
they consider the dual Hamiltonian H{t) for which the 
evolution operator is U'(i). Evidently 

H(t) = -U\t)i?-U{t) = -U\t)h{t)U{t). (2) 



Tong et al. now formulate the putative inconsistency 
as follows: First they argue that H{t) satisfies the con- 
ditions for the adiabatic theorem as well as does hit). 
Then they demonstrate that the adiabatic approxima- 
tion can be very inaccurate for the dual Hamiltonian. 
These conflicting observations constitute the "inconsis- 
tency" . Marzlin and Sanders Q originally formulated 
the inconsistency in a different form; below we will 
show the equivalence of the two formulations. 

In this Letter we examine the dual Hamiltonian H (t) 
and find it violates well known adiabaticity conditions; 
hence there is no inconsistency. The chief difficulty 
in determining whether H(t) satisfies adiabatic condi- 
tions is that we do not have an explicit expression for 
H{t) except in the special cases where the dynamics 
of h{t) are simple enough to allow evaluation of U(t). 
Nonetheless, we are able to give a general argument 
that Hit) violates the conditions of the adiabatic the- 
orem. As an illustration, we apply our general argu- 
ments to a solvable two-level model also studied by ref 
0. In this case it is possible to obtain the explicit 
form of Hit); inspection of this form is sufficient to 
show immediately that H (t) violates the conditions of 
the adiabatic theorem, consistent with our result. 

It is helpful to first review the conditions under 
which the adiabatic approximation is accurate. Con- 
sider a two-level system with the Schrodinger equation 

d ( ci A _ 1 / ujo ne- 1 ^ \f ci \ , . 
dt{c 2 J- 2\ne^ -wo ){c 2 J- W 

We regard the off-diagonal terms as a perturbation 
and ask when the perturbation may be neglected. 
Evidently, the perturbation must be small in magni- 
tude, but even if it is, it can have a big effect on 
resonance when lj w loq. This is more transparent 
if we go over to the interaction picture by writing 
ci = ai exp(— iu>ot/2), c 2 — a 2 exp(iwoi/2). In the in- 
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Away from resonance the off-diagonal terms oscillate 
rapidly and average to zero. Near resonance, the 
perturbation varies slowly in the interaction picture 
and can have a big effect, producing Rabi oscilla- 
tions. The precise condition to be off-resonance is 
\u — ujq\ 3> f2. For this problem, the condition can be 
derived by transforming to another rotating frame via 
a\ = bi cxp[— i(u> — uJn)t/2], a2 — 02 exp[i(w — u>o)t/2\. 
In this frame, the Hamiltonian is time independent and 
equal to (1/2) [(w — ujq)o- z + £Io~ x ]. In summary, we need 
the perturbation to be off- resonance (Q <C \lu — luq\) for 
it to be truly negligible. 

Now let us consider the problem of a general time de- 
pendent hamiltonian h{t). It is helpful to use a slightly 
different approach from the one above, which is useful 
for adiabatic perturbations. To this end, we introduce 
the instantaneous eigenstates \n{t)) that satisfy 



h{t)\n{t)) = e n {t)\n{t)). 



(5) 



We choose the phases of \n(t)) to satisfy 
(n{t)\d / dt\n{t)) — 0, a convention called the parallel 
transport gauge We expand the state of the 

system in this time-dependent basis. Thus 
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In this moving frame the time-dependent Schrodinger 
equation has the form 



i^ t <t> n {t) = ^ A nm (t)exp(z J dt'[e n (t')-E m (t')])(t> n (t). 

(7) 



Here 



A nm (t) = -i(n\ — \m). 

at 



(8) 



The adiabatic approximation amounts to neglect of the 
terms on the right hand side of eq 10. For this to be 
justified, by analogy to eq we see that the neglected 
terms must be off-resonance. Roughly the condition to 
be off-resonance is that the neglected terms should not 
vary slowly; i.e. the terms should be small in magni- 
tude compared to their predominant frequency. Note 
that this is somewhat different from the adiabatic con- 
dition \A nm \ <C \s n — e m \ used by refs 0, @]. More 
precisely, if the Hamiltonian varies predominantly at a 
frequency uj, and the typical spacing between instanta- 
neous eigenvalues is A, our adiabatic condition states 



|^4nm| <C (A — uj). Alternatively, if the Hamiltonian 
varies on a time scale T for which 1/AT -C 1, and 
A nm ~ 1/T, the correction to the adiabatic approxi- 
mation is ~ 1/AT which is very small. This condition 
is more restrictive than neccessary, but it agrees with 
the intuitive expectation that for the adiabatic approx- 
imation to apply, the Hamiltonian must " vary slowly" . 

Sampling the literature we find that the graduate 
text by Schiff [12j gives a rather complete discussion 
of the adiabatic approximation, emphasizing that the 
neglected terms must be non- resonant. On the other 
hand, Landau and Lifshitz |l3j give the more restric- 
tive adiabatic condition T — > 00. Moody, Shapere and 
Wilczek |14| compute non-perturbative corrections to 
adiabatic evolution using 1/AT as the small param- 
eter. Berry 0] also regards 1/AT as the adiabatic 
parameter and suggests that the corrections vanish as 
exp(— AT), also proposed by Hwang and Pechukas (lif . 
Thus it appears to be generally accepted that it is suf- 
ficient for the Hamiltonian to be slowly varying, but 
less restrictive conditions are also discussed. 

For the dual Hamiltonian H(t), denote the instanta- 
neous eigenstates \n(t);H) and the eigenvalues £%(t). 
With the parallel transport gauge, in the adiabatic 
frame 



\m) = 



£^(*)expH f 'dt'E%(t')}\n(ty,H). (9) 



By analogy to eq (JZJ) the Schrodinger equation obeyed 
by the amplitudes <f>^(t) is 



< ro (*)exp(i fdt'[e« (t')-eg(t')])d>»(t). 

_L JO 
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Although, as noted above, we do not have an explicit 
expression for H(t), it is easy to relate the eigenstates 
of H to those of h. Evidently W(t)\n(t)} is an eigen- 
state of H(t) with eigenvalue —e n (t). Hence we write 



(11) 



i(t);H) = C/ t (t)|n(t))expH / dt's n (t% 



The phase of \n(t)\H) has been chosen to ensure par- 
allel transport, (n; H\dt\n; H) = 0. To verify this, it is 
helpful to recall that illdtW — —h, which follows from 
the Schrodinger equation (JTJ. Using eq and this 
result, one can show that 

A% m (t) = A nm (t)ex.p(i[ dt'[e m (t')-e n (t')}). (12) 



Substituting eq l|12|) into eq (|10fl we obtain the final 
form of the Schrodinger equation in the adiabatic frame 
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for the dual Hamiltonian H(t), 



i M^(t) = E A nm {t)<t>l{t). (13) 



Eq (|13|l is the main result of our general analysis. 

Eq 1|13|) shows that the terms that would be ne- 
glected in the adiabatic approximation A nm (t) vary 
slowly. The typical frequency of these terms is 1/T, 
the same as their magnitude. Thus H (t) does not fulfil 
the adiabatic condition, and using the adiabatic ap- 
proximation for H(t) leads to inaccurate results. 

This concludes our general analysis of H(t). We turn 
to a solvable example. Take the two-level Hamiltonian 



cos (9 sin6>e" 



(14) 



Physically we can picture this as a spin i particle in 
a magnetic field tilted at an angle 9 to the z-axis and 
rotating at a frequency to. Essentially this model was 
also studied in ref || • The instantaneous eigenvalues of 
this Hamiltonian are e±(t) = ±wo/2; the correspond- 
ing eigenspinors in the parallel transport gauge are 



! + (*)> = ( 21%t )exp[-iy(l-<-«.s«)]. 



I - (*)> = 



9_ p — iujt 

1 exp[i — (l-cos0)](15) 



A straightforward computation reveals that 

A+_(t) = (+\dtH = |sin0e— tcos9 . (16) 

It follows from eq J7| that the adiabatic frame 
Schrodinger equation for h(t) is 



d 



A* + _e 



A+^e iuJot 

iu t Q 



(17) 

Clearly for the off-diagonal term to be off-resonance, we 
need |o;cos^ — Wo| 3> uism.9. Evidently, this condition 
is satisfied if u <C ljq, i.e., hit) varies slowly. 

It follows from eq (|13H that the adiabatic frame 
Schrodinger equation for H(t) is 







o 
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(18) 



Clearly, for the off-diagonal term to be off-resonance 
we need |wsin#| <C |wcos0|, a condition met only when 
0<1 (or7T-0<l). 

In summary we find that for h(t) to be adiabatic it is 
sufficient that w«wo; but for H(t) to be adiabatic we 



also need 9 <C 1. In this limit, 9^0, the fidelity of the 
adiabatic solution to H(t) (the overlap of the adiabatic 
and exact solutions) computed by 0, approaches 
unity, consistent with our finding. 

We note that in this solvable problem it is possible 
to compute U and explicitly obtain H(t). This column 
is far too narrow to write the entire expression, but it 
includes terms that oscillate at a frequency 



U)q + w 2 + 2luquj cos ( 



(19) 



In the limit wCwo, needed for the Hamiltonian h(t) to 
be adiabatically varying, v — > ujq. Thus, even without 
going to to the adiabatic frame, a cursory inspection of 
H{t) is sufficient to show it is not slowly varying and 
is unlikely to satisfy the adiabatic condition. 

We briefly comment on cases where the adiabatic 
approximation applies to h(t) and H(t). This happens 
when there is a parameter other than T which can 
be tuned to make the magnitude of the off-diagonal 
term in eq (| 1 31) small compared to the predominant 
frequency 1/T. However, for h(t) the approximation 
becomes more accurate as T — * oo (with exponen- 
tially small corrections according to refs [lil llil flfij|). 
whereas for for H{t) the distance from resonance is 
essentially independent of T and is controlled by the 
additional parameter. This can be seen in the exam- 
ple above, where T — > 2ir/oj and 9 is the additional 
parameter. 

We make a few observations in passing here about 
the discussion of the adiabatic conditions after Eq.JSJ). 
First, the effects of A nm are implicitly integrated over 
a finite time window. This is appropriate, since in 
experiments h(t) is typically varied only within a fi- 
nite time window. However, it is the discontinuity in 
dh/dt that causes the ~ 1/(AT) correction to the adi- 
abatic approximation for large T; if all derivatives of 
h(t) are continuous for all t, it is known 0,0,0] that 
the correction will be exponentially small in TA. Sec- 
ond, A was effectively taken to be time independent. 
When both A nm and A vary sinusoidally at frequency 
uj, higher order resonances result when u> = Ao/fc for 
integer k. This is essentially the same mechanism that 
causes resonances going beyond first order perturba- 
tion theory. 

Finally we discuss the equivalence between the for- 
mulations of refs and ||. They assume that the 
adiabatic approximation is accurate for h(t) : U(t) w 
E4dia(i) where 



EW*) =5^|n(t))(n(0)|raq> 



-i / dt'e n {t') 



(20) 



is the adiabatic approximation to the exact evolution 
U(t). Marzlin et al. then develop an approx- 
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imation to U'(t) that we denote V f (t). They ap- 
proximate the dual Hamiltonian [defined by eq (0] as 
#adL(*) = -Ul dia {t)h{t)U &d Ut). They compute the 
evolution V\t) corresponding to H a ^ ia (t), obtaining 



( 2) 

and H adia differ by a small resonant perturbation lead- 



V rt W=^|n(0))(n(0)|exp 



* / dt's n {t') 



(21) 



In contrast, Tong et al. compute a different approx- 
imation to U'(t) that we denote W'(t). They work 
with the exact dual Hamiltonian H(t), but work out 
its evolution using the adiabatic approximation. Anal- 
ogy to eq (|2U[I and use of eq leads to the result 



WHt)=Y,UHt)\n(t))(n(0)\. 



(22) 



The operators and W* are different in appearance 
and in the approximations that lead to them, but we 
will show they are equivalent to the extent that the 
adiabatic approximation U ~ J7 a dia is valid. 

With this notation established we now turn to the 
inconsistencies. Marzlin and Sanders Q consider 
the identity UW = 1 and replace U — + C4dia and 
W — * V* with the disastrous result that UW — ► 
J2 n \ n (t)) ( n (®)\ Tong et al. derive the same in- 

consistency by replacing W with W' : UW — > UW^ 
and using Ea. (|22|l . The equivalence of the two ap- 
proaches can be seen by starting with the trivial iden- 
tity = U^UU^ . If on the right hand side we replace 
WUW -» WUzdizV* 1 we obtain by use of eqs fflty. 
(I21|l and l|22|l . In other words, the approximation of 
Marzlin and Sanders with the adiabatic approximation 
is equivalent to the approximation of Tong et al. 

The resolution of the inconsistencies of Refs.0,0 is 
that ~ W and ~ are not good approx- 
imations; the adiabatic approximation J7 a< ai a ~ U is 
not at fault. In this Letter we have explained the fail- 
ure of the approximation W' ss U* . This may also 
be considered a resolution of the Marzlin and Sanders 
form of the inconsistency, due to the equivalence shown 
above. An alternative resolution was provided in ref 
who used the adiabatic approximation to calcu- 
late a second approximation to the dual Hamiltonian 

-^adiaW = -^adia^^adia- % construction, the evo- 



lution operator corresponding to H adia (t) is iX[ dia (t) 
which is indeed a good approximation to U* . Since the 
trustworthy adiabatic approximation has been used to 
compute both H adia (t) and H^ ia (t) we may presume 

that H^Jt) « H^M « H(t). This does not mean 
that the three Hamiltonians will generate essentially 
the same evolution since small errors can grow upon 
time exponentiation. Indeed ref Q finds that H ad \ a 



ing to the inference that may be very different from 
^adia- Since Ul dia rj W this may be considered to ex- 
plain the failure of the approximation V' w UK 

To summarise, we have studied the dual pair of 
Hamiltonians h(t) and H(t) that generate time evolu- 
tion U{t) and U'(t) respectively. Marzlin and Sanders 
showed that if an adiabatic approximation is used 
to compute H{t) an inconsistency results. Tong et al. 
showed that the same inconsistency results if the 
adiabatic approximation is applied to the exact dual 
Hamiltonian H(t). We show that these observations 
are essentially equivalent. Our main finding is that 
even if h(t) satisfies the conditions of the adiabatic 
theorem, H(t) will not (except in the trivial case when 
the total change in h(t) is small) because the terms 
neglected in the adiabatic approximation are resonant 
for H(t). Thus the breakdown of the adiabatic approx- 
imation for H is not inconsistent with the adiabatic 
theorem. 

It is a pleasure to acknowledge helpful discussions 
with Francesc Ferrer. 
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